Robust optical delay lines via topological protection 
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Phenomena associated with topological properties of physical systems are naturally robust against 
perturbations. This robustness is exemplified by quantized conductance and edge state transport in 
the quantum Hall and quantum spin Hall effects. Here we show how exploiting topological properties 
of optical systems can be used to implement robust photonic devices. We demonstrate how quantum 
spin Hall Hamiltonians can be created with linear optical elements using a network of coupled 
resonator optical waveguides (CROW) in two dimensions. We find that key features of quantum 
Hall systems, including the characteristic Hofstadter butterfly and robust edge state transport, can 
\ be obtained in such systems. As a specific application, we show that the topological protection can 



o 



be used to dramatically improve the performance of optical delay lines and to overcome limitations 
, related to disorder in photonic technologies. 



Particles in two-dimensional structures with a magnetic field exhibit a remarkable variety of macroscopic quantum 
phenomena, including integer [l| and fractional EJ quantum Hall and quantum spin Hall effects and predicted 
regimes of fractional or non-abelian statistics 0, [J A hallmark of these systems is the presence of edge states, whose 
O transport properties are robust against disorder and scattering, leading to quantized conductance sufficient to provide 
+L \ a resistance standard 0, Q] ■ Natural robustness of topological states is actively explored in quantum computation 
0, Q . Recently, approaches to observin g si milar quantum Hall behavior in bosonic systems including ultra-cold gases 
^ . (for a review see Ref . fiof ) and photons [llT - [lti | have been suggested. 

Our method for realization of topological protected photonic devices makes use of two dimensional arrays of coupled 
resonator optical waveguides (CROW) to simulate a 2D magnetic tight-binding Hamiltonian with degenerate clockwise 
and counter-clockwise modes. This approach does not require explicit time-reversal symmetry breaking 11 -0, but 
\^ , the degenerate modes — time-reversed pairs — behave analogously to spins with spin-orbit coupling in the electronic 
quantum spin Hall effect (QSHE) 17T4l9jj . and experience a spin-dependent magnetic field (Fig.Q}. When the clockwise 
and counter-clockwise modes are decoupled, we can selectively drive each mode and observe quantum Hall behaviors 
without breaking the time-reversal symmetry in the tight binding Hamiltonian. In a direct analogy to the electronic 
integer quantum Hall effect, we show that photonic edge states carry light at the perimeter of the system, while 
being insensitive to disorder, and therefore forms a basis for robust photonic devices. In particular, in comparison 
to state-of-the-art 1-D CROW systems, our approach can be dramatically more resistant to scattering disorders and 
^ , fabrication errors. 

X: 
ft ■ 

2D Photonic System and Quantum Spin Hall Hamiltonian 



As illustrated in Fig. [TJ our system comprises optical ring microresonators that support degenerate clockwise 
and counter-clockwise modes, restricted to one pair per resonator. We consider these modes as two components 
of a pseudo-spin, i.e., clockwise (a — —1, or psuedo-spin down) and counter-clockwise (er = +1, pseudo-spin up) 
circulation. Resonators are evanescently coupled to each other and have been studied in the context of ID CROW 
[2f3 |. where the coupling leads to a tight-binding model for photons and the corresponding photonic band structure. By 
coupling these modes in a two-dimensional arrangement, as we show below under appropriate conditions, the dynamics 
of such photonic system is described by a Hamiltonian for charged bosons on a square lattice (tight-binding), but 
with the addition of a perpendicular, spin-dependent effective magnetic field: 
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Figure 1: Schematics of the photonic system with a synthetic magnetic field: (a) The dynamics of two coupled 
resonators is governed by a magnetic-type Hamiltonian (equation 3). The lengths of the upper and lower branch differ from 
each other so that the phase difference is 4na and m is an integer, (b) A magnetic tight-binding model can be implemented in 
a 2D lattice of coupled resonators. Only the waveguide phase differences are shown, (c) By introducing scatterer in resonators 
and waveguides, one can induce in-plane magnetic field (Zeeman) and spin-flip hopping terms, in analogy to quantum spin Hall 
models. 
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where k is the tunneling rate of optical modes and a\ x y is the photon creation operator at resonator at site (x,y) 
with different spin componenets (a = ±1). Specifically, photons acquire a 2iraa phase when they go around a 



plaquette-equivalent to having a quanta of magnetic flux penetrating each unit plaquette 2l|, |22j . 



First, we show how the description of our system using such a Hamiltonian is derived. We start by considering two 
coupled resonators (Fig. [TJa)), focusing only on the counter-clockwise modes inside each resonator. The length of 
connecting waveguides is chosen such that photons destructively (constructively) interfere inside the waveguide loop 
(resonator), respectively, and therefore, they will be confined in the resonators rather than waveguides. Moreover, the 
lengths of the upper and lower branches of the waveguide differ from each other, so when a photon hops from the left 
to the right resonator, it acquires a different phase than when it hops in the opposite direction. This can be formally 
verified using the standard input-output formalism [23I ]. In particular, the boundary condition at the left resonator 
can be written as 



(2) 



3 



where E x 's are waveguide field operators at the vicinity of the x-th resonator and a x is the resonator field operator, as 
shown in Fig. HJa). The resonator field equation of motion is: d t d x — —na x — ^/2~K,E x n , and similarly for the right res- 
onator. Photons propagate freely between the resonators, so for the upper branch we have E%\ x = —iE° ut exp(— 2iria), 
and similarly for the lower branch. By eliminating the waveguide fields, the left and right resonator fields dynamics 
will be given by $ta x ( x _|_i) = iftexp (±27ria)a x+ i( x ), consistent with tunneling between resonators. Therefore, the 
effective Hamiltonian of the two resonators takes the form: 

Htwo-res = -Ka t x+1 a x e- 27 ' la - Kala x+1 e 27 ' la . (3) 

The above analysis for counter-clockwise modes (pseudo-spin up a,f x ,y) in the resonators shows that, in the absence 
of backscattering, they are decoupled from their time-reversed counterpart, i.e., the clockwise mode of the resonator 
(pseudo-spin down aj, X)J/ ). At the same time, the pseudo-spin down component will experience a magnetic field similar 
to the pseudo-spin up component, where only the sign of magnetic field is changed (a —a). Now, by connecting 
resonators in a lattice structure and tuning the phase of the connecting waveguides, we can arrange the acquired phase 
around each plaquette to be uniform and equal to 2ira. The phase can be tuned either by changing the length (or the 
index of refraction) of the connecting waveguides or by coupling ring resonators to the sides of the waveguides (similar 
to a Mach-Zender configuration 24j,|25|). The implementation of a Landau-type gauge is shown in Fig.QTb) where the 
corresponding Hamiltonian is the form of equation ((T|). Indeed, this is the trivial form of Hamiltonians with QSHE in 
an analogy to time-reversal invariant spin-orbit interactions in solid state systems 17-1^], where the magnetic field 



has opposite orientations for the two pseudo-spin components. Furthermore, one can control the coupling between 
different pseudo-spin component and exploit a wider class of QSHE Hamiltonians on a square lattice. In particular, 
semi-transparent scatterers inside the resonators or the connecting waveguides can be engineered to mix different 
pseudo-spin components with each other. 



Generalized pseudo-spin-orbit interaction 

To illustrate this mixing, we consider the addition of a pair of scatterers in every vertical connecting waveguide 
as shown in Fig. [lie). For simplicity, we assume the scatterer is weak and loss- less. We characterize the strength of 
the scatterer by a parameter e, where the transmission coefficient is near unity (t s ~ 1) and the reflection coefficient 
is r s = ie/y/2. As shown in the Supplementary Information (SI), the corresponding Hamiltonian of a single vertical 
array of resonators will be: 



u m = -*E( 4™ +1 aLv+i ) T v + h.c. (4) 




The diagonal terms are identical to the tight-biding terms of equation (JTJ and the off-diagonal terms represent 
hopping between two adjacent sites while undergoing a spin- flip, i.e., the scatterers couple clockwise to counter- 
clockwise photons (Fig. [He)). This spin-flip hopping is similar to the Rashba term in the context of spin-orbit 
interaction 
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Similarly, if we consider a pair of weak scatterers inside the resonators, then corresponding Hamiltonian takes the 
form: 




+ h.c. 



(5) 
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The first term is the usual tight-binding form and the second term represents the in-plane magnetic field which is 
enhanced by the finesse of the resonators (i.e., number of photon round trips T ~ 7r/(l — f 2 ))- If these vertical arrays 
replace the vertical arrays of Fig. [ljb), then the overall Hamiltonian of the system encompasses both an in-plane 
Zeeman term, perpendicular to the synthetic magnetic field (due to on-site scatters) and a hopping-spin flip term 
(similar to Rashba interaction). 



Probing the photonic system 

We now show how optical spectroscopy measurements can be harnessed to analyze the transport properties of our 
photonic system. In particular, as shown in Fig. [11(b), by coupling two waveguides to the lattice edges, we evaluate 
the transmission and the reflection of an input light field and study their properties in connection with the magnetic 
states of the system (localized states and edge states). 

We first consider a trivial case of quantum spin Hall effect where pseudo-spin-flip terms are absent. In this regime, 
we have two decoupled copies of regular quantum Hall states for opposite pseudo-spin components. We restrict our 
analysis to a single spin component, and for brevity, we drop the spin index. This choice enables us to examine in 
detail methods for probing the system and determining its response to errors and disorders, without the additional 
complexity of spin Hall physics. 

Using a formalism similar to the quantum scattering theory, we investigate the problem of scattering of light field in 
optical waveguides connected to our photonic system and evaluate transmission and reflection coefficient under various 
conditions. The waveguides only couple to co-propagating modes in the resonators (counter-clockwise in Fig. Q3b)), 
and thus under our assumption, the reflection in the input channel and also transmission in the output channel is 
zero (i.e., R,T' = shown in Fig. [11(b)). The input-output probing waveguides are coupled to two resonators in the 
systems denoted by \in) and \out), respectively. As shown in SI, the self-energy of these resonators can be written as 
S = — i%\in)(in\ — i%\out)(out\, where the coupling strength is defined as v. Using Lippman-Schwinger equation, one 



can deduce different reflection/transmission coefficients 271 l28j . In particular, the reflection coefficient is given by 



r'(w) = -iv 

zEin,j€zout 



H 



(6) 



which means an appreciable reflection should be observed when the frequency of an incoming photon becomes resonant 
with the energy of a photonic state inside the system. Note that if the photonic system is a single resonator, equation 
(6) reduces to the familiar form of: r'(to) = — 7^ — r. 

The energy spectrum of the Hamiltonian of equation (JT]) for an infinite lattice is the well-known Hofstadter butterfly 
0. We consider a N x x N y lattice with torus boundary condition (i.e., coupling top-bottom and left- right edge 
together) to simulate the effect of an infinite lattice. According to Hofstadter spectrum, for rational magnetic fluxes 
(a = p/q), each magnetic band has many states ( Nx y -fold), which is the reminiscent of the Landau degeneracy in the 
continuum. The result of our numerical solution is shown in Fig.[2]where the reflectivity (i?' = |r'(w)| 2 ) is evaluated 
for different frequencies and magnetic field (a), by the formalism described above. High reflectivity occurs when the 
lower waveguide light is coupled to the system and completely transferred to the reflection output channel (the second 
waveguide), similar to a channel drop filter. We can readily see that the energy spectrum of the uncoupled system 
(Hofstadter butterfly) can be obtained by measuring the system reflectivity. We note that in order to to resolve 
different energy levels in the spectrum, v should be chosen to be sufficiently narrow (< 8k/(N x N v )). 



Photonic edge states 



In a direct analogy to 2D electrons in a magnetic field (quantum Hall physics) 29- 31), we recognize quasi-one 



dimensional states that are localized at the perimeter of the system which carry current and are immune to disorders 
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Figure 2: Hofstadter butterfly spectrum: Each point represents a reflectivity greater than 0.005, for a 10x10 lattice with 
torus boundary condition and the coupling i//k= .02. The green line shows the spectrum at the magnetic field of interest (a) 
for the rest of the figures in the article. 



in form of random potential. In contrast to toroidal boundary conditions where only magnetic bands are present, in a 
finite square lattice, there exist states between magnetic bands which are known as "edge states" . These edge states 
carry a chiral current around the system perimeter, as shown in Fig. G^a). For each edge state, there is a similar edge 
state with an opposite chirality which has the same energy magnitude but with the opposite sign. 

So far we have considered a perfect system. To illustrate the robustness of the system to disorder, we consider 
a resonator that is detuned from its neighbors provides a model for "non-magnetic" disorder. This disorder can be 
readily characterized by an on-site potential at that site {Ua\. y a Xyy ). The main source of such frequency mismatch is 
small variations of the waveguide widths and the perimeters of the microrings during the fabrication process. Such 
imperfections are a common problem in photonics and prevents coupling large number of resonators 32|, [33 1 . 

In electronic quantum Hall systems, the edge states are immune to disorder Q. We find that such robustness 
applies to our photonics system. In particular, when such a disorder is located 'inside' the edge state, the edge state is 
obviously not affected. However, when the disorder is located on the edge, the edge state routes around it, as shown 
in Fig. [3Jb) . More precisely, scattering which would reverse the current is prevented due to the lack of additional 
states with appropriate energy and momentum. 



Application to delay lines 

The transport properties of edge states provide a robust alternative to conventional CROW in photonic delay lines. 
In particular, we compare the transport properties of our photonic quantum Hall system to CROW, as illustrated in 
Fig. 2J In the quantum Hall system, there is a robust transfer band provided by edge states which carry photons from 
the input waveguide to the output waveguide (Fig. [4] inset), in a direct analogy to electronic edge state in the context 



of the integer quantum Hall effect [29 1. 



In both systems, the operational bandwidth is given by smooth, linear part of the dispersion relation. In the 
quantum Hall system, the edge state band (Fig. [3]), is located between two Hofstadter bands (Fig. [5), while in the 
CROW configuration, the operational bandwidth is in the middle of the tight-binding dispersion to avoid the group 
velocity dispersion [20j. Moreover, in both systems, the delay time is proportional to the number of resonators involved 
in the transport. In the quantum Hall system, the transport can be either performed along the long or short edge 
of the system, depending on the input frequency, as shown schematically in Fig. [4] inset. Therefore, in both systems 
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Figure 3: Edge states and their dispersion: Probability current for an edge state in the absence (a) and in the presence 
of a disorder (b). The location of the disorder is shown by the black dot, and U/k = 5. (c) shows the dispersion relation, a 
wave number is evaluated for each energy eigenstate at the lower (red) and upper edge (blue), i.e., KA is the phase difference 
between two consecutive resonator at the edges. Magnetic bands state do not have well-defined wavenumber and are shown in 
the yellow area. The dashed line shows the state corresponding to (a). In these plots, a 10x10 square lattice with a = \ is 
considered. 



in the absence of disorder, the bandwidth-delay product increases by the length (perimeter) of the system. However, 
in the presence of disorder, CROW and edge state behave differently as the system size increases. In particular, in 



ID systems (e.g. CROW) the disorder leads to localization 34j-436jj and therefore, the transmission is impeded. In 
contrast, the transport of edge states is topologically protected against disorder 0, |3t| and the transmission is not 
affected (as discussed above for Fig. [3]). 

To confirm that our edge states provide a robust transport, we numerically study the effect of disorder, by assigning 
a random frequency mismatch to the resonators. Taking average over many different disorder realizations, we observe 
that by increasing the size of the system and consequently the bandwidth-delay product, the transmission in CROW 
decreases while the transmission with edge states is unaffected, as shown in Fig. [5l Furthermore, we observe that 
while magnetic band states and CROW depend sensitively on the disorder (position/strength), from one realization 
to another, the edge states are insensitive to the specific parameters, as shown in the standard deviation of the 
reflectivity and delay time in Fig. 0Ja,b). The effect of the disorder on edge states is simply a shift of energies, which 
is manifested in a frequency shift in the edge states transfer band. 

We note that as time-reversal symmetry is not broken in our system, we can not utilize such edge states as a one-way 



waveguide similar to Ref.[38(. More precisely, when we input a light field in the backward direction into the system 
(by swapping the input and output channel), the waveguides couple to the opposite rotating field in the resonators 
(opposite pseudo-spin) and experience a magnetic field with an opposite sign. Therefore, the system is reciprocal and 
the transport properties of the forward and backward feed are identical to each other. 

Finally, we investigate the effect of loss in the resonators and other imperfections. Loss can be represented by a 
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Figure 4: Edge states vs. CROW: This figure compares the quantum Hall system (left panels, 10x10 lattice) and CROW 
(right panels, array of 40) performances as delay lines. Blue curves show the average time delay (output reflectivity) in a(b), 
respectively, and the gray area highlights the standard deviations in the presence of non-magnetic disorder (a gaussian disorder 
with a width U/k = 0.4 for 500 realizations). While the transport is quite noisy in the magnetic bands and CROW (tight- 
binding band), the edge state bands exhibit noiseless transport with delays comparable to CROW. Depending on the input 
frequency, different edge states participate in the transport which leads to shorter or longer delays, as shown in the insets. In 
the presence of magnetic disorder (the strength is a gaussian with a width ej- = 0.1 and the phase is random [0, 2tt)) and loss 
(/tin = 0.02k), the transport properties are more degraded for CROW than edge state bands. While the counter-clockwise 
modes of the resonators is excited through an input field, the onsite scatterers backscatter photons in the clockwise modes. 
These modes leak out into R and T' channels which are non-zero in these plots. The coupling between input-output waveguides 
and the system is chosen to optimize the transport (for edge states v = 6k and for CROW v — 2k). In the quantum Hall 
system, the input and output waveguides are coupled to (x = 2, y — 1) and (x = N x — 1, y = 1) resonators, respectively. 



non-hermitian term in the Hamiltonian: —ini n a\.a x where the dissipation rate is given as the product of loss per 
roundtrip, number of roundtrips per decay time and the photon decay rate [K% n = /u(^-)(4k)) (see SI). The photonic 
loss can attenuate the reflection in the edge state transfer band due to the propagation around perimeter (Fig. S£c)). 
Silicon-on-insulator technology, where up to 100 micro-ring resonators have been successfully coupled to each other 
33], is a promising candidate for implementation of our scheme (see SI for experimental relevant parameters). 

Other types of imperfection such as surface roughness can cause undesired backscattering which mixes pseudo-spin 
up and down, acting as "magnetic disorder" . As show in the SI, these imperfections can be modeled by a magnetic 
disorder Hamiltonian. The backscattering effect manifests in the reduction of signal in R',T channels and some 
leakage in R, T' channels. Fig. |4jc) shows different transport coefficients. We observe that although the transport 
properties of magnetic bands states are affected by such magnetic disorder, edge state transport is robust, which is 
due the suppression of backscattering events. In particular, the scattering of a forward-going spin-up into forward 
(backward) -going spin-down is inhibited due to energy (momentum) mismatch. 
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Figure 5: Transmission in ID CROW and 2D quantum Hall system: For ID CROW (2D quantum Hall), the transmis- 
sion (R') is evaluated at the center frequency uj — (u — 1.5k), where the group velocity dispersion is minimum, respectively. 
A gaussian disorder is assumed with a width U/k = 0.4 and the transmission is averaged over an ensemble of 500 realizations. 
Gray shades show the standard deviation. By increasing the bandwidth-delay product (x-axis), the transmission for ID CROW 
decreases and the noise increases, while both the transmission and the noise in the 2D quantum Hall system remains constant. 



Outlook 



Optical signals might be a promising alternative to electronic signals in the future circuits [391 ] . One of the key 
requirements is the ability to filter and slow down light on-chip over a large bandwidth (several Gbps) for various time- 
domain processing such as optical buffering and multiplexing 4CJ. However, the effect of disorder in the millimeter 
size footprint is detrimental, e.g. unwanted signal modulation in transmission spectrum 33]. Our system provides a 
platform to realize photonic system immune to disorder. Potentially, one can envisage using other types of topological 
properties to improve photonic technologies. 

In addition, our photonic system enables a new approach for exploration of various fundamental quantum Hall 
phenomena. This photonic system not only enables investigations of quantum Hall physics by simulating different 
types of Hamiltonians at room temperature, but it also taps into topological features to provide robust devices for 
photonics. In the non- interacting regime (which was the topic of this article), one can explore the Hofstadter butterfly 
of photons, and photonic edge states as delay lines immune to disorders and also localization in 2D for non-interacting 
particles 



41 



Furthermore, with the addition of interaction between photons, this system opens up exciting prospects 
for exploring many-body topological state of light. Although the ground state properties of such systems have been 
extensively studied, the suitable characterization and measurement of strongly interacting photons is still an open 
question. In particular, the chemical energy is absent for photons and the relevant conditions to study photons 
involves an externally driven system which leads to non-equilibrium situations 42[. Another advantage of photons 
is their flexibility to form various system topologies (torus with different genera) by simply connecting waveguides 
to each other and manipulate such states for topological quantum computation. Intriguing additional applications of 
these ideas await further exploration. 
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Supplementary Information 



I. TRANSFER MATRIX FORMALISM AND TIGHT-BINDING MODEL 



Here we derive the magnetic tight-biding model which is presented in the main text as equation (1). In contrast 
to the toy model of two coupled resonators, presented in the main text, we consider an infinite array of coupled 
resonators similar to CROW [2(3]. 0m system consists of an array of coupled resonator and waveguides which is 
periodic and linear. Therefore, the transfer matrix formalism and Bloch's theorem are suitable to treat the system 
dynamics 20l. l43l. l44j . in contrast to the input-output formalism presented in the main text. 

We use the transfer matrix formalism for a ID array of coupled resonators. As shown in FigJH the fields at the 
right edge of the n-th resonator and the n-th waveguide are respectively characterized by: 



f — 



( a n \ 




(<\ 


K 


£' - 


b' n 












\ d 'nJ 



(7) 



such that the electric field inside the resonators is: 



E n {p,4>) = /(/>)* 



a n exp 



6„exp 



ij3R( 



7T — (j)^ 



-i/3R{ 



d n exp 



-i/3R{ 



77 ~ 4>\ 







+ c n exp 






7T 



for < 4> < k 
for -7r < 6 < 



(8) 
(9) 



where p is the radial coordinate, <j) is the azimuthal angle relative to the counter-clockwise propagation direction and 
f(p) is the radial dependence of the electric field. 

We assume the resonators have length 2R and the total length of the connecting waveguides is 2L, as shown in 
FigEl Therefore, transfer matrices corresponding to free propagation inside resonators and waveguides are: 




Figure 6: ID array of coupled resonators and waveguides 
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0R 



M res (R) 



M wg (L) = 





-i/3R 





e i/3R 
\ e 

/ e i/3L+i2TTay q 






-i/3R 



\ 



-iftL-{-i27ray 

o 
o 






^ij3L—i2-Ka.y 

o 



(10) 



(ii) 



respectively, where the wave number is j3 — lj/c and uj is the light frequency. We assume that the clockwise and 
counter-clockwise modes, i.e., forward and backward fields, are decoupled from each other, which sets the off-diagonal 
elements in M wg and M res to be zero. The regime where this coupling between forward and backward fields becomes 
important is discussed in the main text and also in Sec lIVI of this supplemental information (SI). We further assume 
that the coupling regions are long enough (compared to the wavelength of light) so that the coupling between clockwise 
and counter-clockwise photons (e.g. the coupling between (a,b) and (c',d') fields) can be ignored. In other words, a 
single coupling region couples only (a, b) to (a',b') and (c,d) to (c',d'). Therefore, a loss-less coupling region can be 
characterized by two parameters t and r, which represent the amount of transmission into and reflection off of the 
adjacent element, respectively, as shown in FigJBl Note that these coefficients obeys the conditions: \t\ 2 + \r\ 2 = 1 
and r*t + rt* = 0. One can easily find the transfer matrix elements of such coupling regions for (a, b) and (a', b'), by 
setting a = 1, a! = £, b = r, b' = and a = 0, a' = r, b = t, V = 1, and similarly for coupling between (c, d) and (c', d'). 
Therefore, the transfer matrix of a single coupling region which couples (a, b) to (a', b') and (c,d) to (c',d') can be 
written as: 



M, 



Cpl 





r 
t 

k 



(12) 



(13) 



In this notation, the field amplitudes in the right edge of the waveguide can be obtained from that of the adjacent 
resonator: E' n — M cp iM res £ n . 

The transfer matrix of the building block of such system {£ n +i — M£ n ) is a product of the above matrices 
M = M cp iM wg M cp iM res . Using Bloch's theorem, the dispersion relation is written as: 

Det[M - exp[iKA)I] = 0, (14) 

where K is the Bloch quasi-momentum and A is the unit spacing. First we solve the case where the phase imbalance 
between waveguide arms is zero (a — 0). We define the detuning from the resonator frequency as A = uj — luq. We 
consider a configuration where resonators are separated by L = (2p + l)Ao/4 where p is a non-zero integer and Ao is 
the resonant wavelength. In this configuration, photon reflections off the resonators interfere destructively with each 
other, and therefore, photons spends most of their time in the resonators rather than the waveguides. We are primarily 
interested in relatively high quality factor resonators. Therefore, we evaluate the above expression to the first order 
in (1 — \r\ 2 ), which is proportional to the inverse of the resonators finesse .F ~ 7r/(l — r 2 )3>l. By assuming r to be 
real and R = n\ for resonators length, Eq. (fT4"l) gives a dispersion relation: cos(K A) = (— l) p ijz$yj~ ^ = (~ 1) P ^' 
where k characterizes the cavity field decay rate from on side of the resonator into the waveguides, i.e., the photon 
decay is 4k = (1 — r 2 )c/R. This decay is the extrinsic desirable decay, in contrast to intrinsic loss (Kin) which we 
discuss in Sec.V. Note that the sign of the tunneling amplitude depends on the parity of p, which we assume to be odd 
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(p = 2m + 1 in Fig. la of the main text), so that we get the conventional minus sign in the tight-binding Hamiltonian. 

Alternatively, one can also diagonalize the transfer matrix of a cell and study its eigenvalues as propagating waves. 
Under the above assumption, the diagonalized form of the transfer matrix of a unit cell will be of the following form: 



M diag (0) = - 



V 



A \2 







'A^2 

. 2k, > 



£-yi-(£) 2 y 



(15) 



We can infer: — ± i^J 1 — (j^) 2 ^ = e ±lkA , which is equivalent to the aforementioned dispersion relation. Again, 
this shows that the propagating wave only exists when |^| < 1. Note that in this situation, we have two modes with 
the same dispersion relation, which corresponds to clockwise and counter-clockwise photons inside the resonators. 

Furthermore, we consider the situation where there is a 27ra phase imbalance between waveguide arms. In this 
case, the diagonalized form of the transfer matrix takes the form: 



M diag (a) = M diag (0) 



V o 

















g — 2Z7TQ 































(16) 



From Eq. (fT4|) one can easily see that the above transfer matrix corresponds to two dispersion relations: cos(-ftf A ± 

_ _ A_ _ u(K)—u 
2k 



2na) = — ^ = — "° ' wmcn is the same dispersion of the tight-binding Hamiltonian in ID with a magnetic field: 



— a crx+l,y a crx,y<i T a 'crx,y a ax+X,y ti 



(17) 



<7,x,y 



In particular, the mode corresponding to (c,d) fields in the resonators (pseudo-spin up) experiences a magnetic field 
(+2na) , while the mode corresponding to (a,d) field (pseudo-spin down) experiences an opposite magnetic field 
(— 2na). In a 2D arrangement, by choosing these imbalances to be a function of rows (±2iray), spin-up (spin-down) 
photons experience a uniform magnetic field corresponding to +27ra(— 2ixa) flux per plaquette, respectively. 



II. DERIVATION OF THE HAMILTONIAN WITH SCATTERERS 



In this section, we show that by implanting a scatterer in the system, one can achieve a photonic version of different 
quantum spin Hall Hamiltonians. In particular, the presence of a scatterer in each resonator leads to an effective 
"in-plane magnetic field" and the presence of scatterer in each connecting waveguides leads to a "spin-flip hopping" 
term. Below, we consider these effects. 



A. Spin-flip hopping 

Here, we study of the effect of scatterers in the connecting waveguides. In particular, we position two scatterers 
symmetrically in the middle of each waveguide, as shown in Fig. 1(c) of the the main text. Therefore, the transfer 
matrix of the scatterer which couples (a,d) and (b,c) field together (as defined in EqJT]), is given by: 
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2 „2 

^ 



i ~ r t 

.2 2 

IT ^ 



M scatt 

ts / 

In this situation, the transfer matrix of the unit cell is given by: 



(18) 



R 



R 



M = M rPX I — ) .M rPS I - ) M cpl .M wg . I - ) .M smtt .M vw [ - | . 



L 



1 scatt- lvl wg 



L 



(19) 



We assume the scatterer is weak so we can write the transmission/reflection coefficients in the following form: 



r s = iy/l - t 2 

t s = l-e 2 / 2 , e < 1 



(20) 
(21) 



In the absence of the scatterer, (a,b) and (c,d) modes couple to each other independently, i.e. (1,0,0,0) couples 
to (0, 1, 0, 0) and (0, 0, 1, 0) couples to (0, 0, 0, 1). As we discussed, these decoupled manifolds form an identical tight- 
biding Hamiltonian for two pseudo-spin components which we call (o-j- , a±). However, the scatterer couples these two 
manifolds to each other. It turns out that a proper basis in the presence of the scatterer is a Hadamard rotated 
basis, i.e. ^(1, 0, ±1, 0) and -^(0, 1, 0, ±1). This can be verified by using a degenerate perturbative method to find 
eigenvalues of the transfer matrix. Therefore, one finds out that the basis which diagonalizes M in the subspace of 
(a,c) or (b,d) is the x-basis -i.e. a^± — -^(a-f- ± aj.). The eigenvalues of the system to the zeroth order in inverse of 
resonator finesse oc (1 — r 2 ) are: 



(2 - ey/T^) ^ ± iy J(2 e 2 ) 2 + (W^ - 4 - Ae 2 + e±){±) 2 
A± - ^ ( 22 ) 

(2 + eV^) ^ ± iy J(2 e 2 ) 2 + (-Ae^A^ - 4 - Ae 2 + e±)(±) 2 
X ± - ^2 • (23) 

These eigenstates have magnitude equal to one. Since each pair of eigenvalues are complex conjugate of each other, 
the four eigenvalues of system can be written in form of Bloch waves: 



A^ = e 
A^ = e 

This leads to two dispersion relations to the first order in e: 



±iKA 
±iK'A 



-(1T^ = cos[XA] 



which can be simply written as: 



uj-ljo = -2/c(l ± e) cos [if A]. 
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Now, we want to write this Hamiltonian in the second quantized form using the original pseudo-spin components 
in z-basis (a^,a^). Note that the above dispersion relation is written in the Hadamard rotated basis (a^,a^) where 
a^± = -^g(a-|. ± oj,). Therefore, we have: 



n i t t \( «(l + e) \ / a_ 



(24) 



The Hamiltonian in the original basis takes the conventional hopping flip term: 




(25) 



<*,j> 



While the direct hopping amplitude is n, the spin-flip hopping amplitude is ne. We note the above Hamiltonian is 
only valid when the scatters are positioned symmetrically in the middle of waveguides. If we move the scatterers 
symmetrically along the waveguides, we arrive at the same form of Hamiltonian with a modified strength of spin-flip 
term. 



B. In-plane magnetic field 

We now evaluate the effect of having a scatterer in every resonator rather than the waveguides. As we show, this 
model results in an in-plane magnetic field term in the Hamiltonian. 

In this case the transfer matrix of the unit cell is: M = M res (R/2)M scatt M res (R/2)M cp iM wg (L/2)M wg (L/2). 
Writing again the transfer matrix in the x-basis, we perturbatively find the eigenvalues to the zeroth order in (1 — r 2 ) 
and the first order in (e): 

cos^A=i(A- + A-) ^ 



2 v + - ' 2k 1 - r 

cos*'A = i(A- + A^-A + _^ 



Therefore, the dispersion relations will be 



4eK 

A = Tt; or - 2Kcos[ifAl. 

(1 - H) 



The corresponding Hamiltonian in the x-basis (a^,a^) is: 



ff =^D«t ( «l, > ( -f + i T )(::)- e < ^ «i, » ( ; :)(: :: ) , ™ 

where the finesse of the resonators T ~ 7r/(l — r 2 ), characterizes the average number of photon round trip before 
leaving the resonator. The above Hamiltonian show that the scatterer inside the resonator lifts degeneracy and splits 
the modes. The splitting is proportional the strength of the scatterer and the number of photon round trips inside 
the resonator. In the original z-basis (a-f,a^), the Hamiltonian takes the form: 
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This means that the presence of the scatterer inside the resonators result in an in-plane magnetic field. The amplitude 
of this in-plane magnetic field is given by the strength of the scatterer and the number of photon round trips inside 
the resonator. Similar to the spin-flip term, the strength of this terms can modified by symmetrically moving the 
scatterers long the resonators. 

In the end, we note that the backscattering can also occur in the waveguides. However, in contrast to the resonator 
case, where the scattering was enhanced by the finesse of the resonators, the backscattering effect in the waveguides 
is not enhanced due the off-resonant nature of the waveguide loops. 



III. INPUT-OUTPUT FORMALISM FOR PROBING THE PHOTONIC SYSTEM 

We now consider the effects of coupling "probe" waveguides to the system. When the system and probe waveguides 
are not coupled to each other the Hamiltonian is given by Hq + H pro i, e . The first term corresponds to the system 
which we want to probe and the the second term corresponds to free propagating light in the probing waveguide, 
i.e., H pro i, e — ^2 k oj(k)\k)(k\ where \k) represents the state of one photon in the waveguide with wave number k. We 
note that since the system is linear and we are interested in the linear response, we treat the system in the first 
quantized representation and do not use creation and annihilation operators. On the other hand, the resonator fields 
are coupled to the forward- and backward-going field of the probing waveguides represented by the following coupling 
to the tight-binding Hamiltonian: 



V = J= J2 Vt\in)(k\+Vt\k)(in\ 

— oo< fc<+oo 

+ -7£ E v^outm+vr'mout] (28) 

— oo<A;<+oo 

\in) (| out)) represents the photonic state of the resonator which is coupled to the input- (or output-) waveguide. 
We assume that resonators to be large enough so that each resonator mode couples to single running wave in the 
probing waveguides. For instance, CCW mode of the resonators is coupled to forward- (backward-) going light field at 
the input (output) probing waveguides, respectively (Fig. 1(b) of the main text). Therefore, only positive (negative) 
momentum terms have non-zero coefficient in input (output) Hamiltonian terms. C is the length of the waveguides 
which enters as a normalization factor. 

Now, we want to study the scattering of an incoming plane wave |fc) which is sent into the probing waveguide. The 
effect of system-probe coupling is manifested in the output light field (\ip)) which leaves the system in the probing 
waveguides. Since the system is linear, one can also look at this process from a single-photon picture: a photon comes 
from the input waveguide, while the system is initially in the vacuum state. Then, the photon interacts with system 
and eventually gets scattered into different output channels. The photon leaves through the probing waveguides and 
the system returns to original vacuum state. 



In order to study such transport effects, we follow the formalism introduced in Refs. 27|, |28j to investigate our 
system. We start by Lippmann-Schwinger equation: 

IV;} - |fc) + * —VM (29) 

^-sys -^probe i &t 



which can be solved through iteration: 
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71=0 



iv>> = £ — t? — ji — --y) i fc >- (30) 

— \Wfc — ri sys — tlprobe + It 1 



In other words, the T-matrix which characterizes the scattering amplitudes into different output channels can be 
written as: 



(k>\T\k) = (k'\"£[ 77 77 --V) \k) (31) 

^ \Wk - Usys - tlprobe + «e / 

Since the total system is linear, each photon should be exchanged an even number of times between the waveguide and 
the system to have a non-zero contribution in the output field. Therefore, only even powers in the above expression 
have non-zero contribution to the T-Matrix. In other words, 

^ 

Tk'k = Sk\k^ ■— -t > Vk> tm G mn {uJk)V n k (32) 

LO k - Wfc/ + it C ' 

where the Green's function is defined as: G — G° sys + G^ygVG^gVG^yg + . . . and G° ys (G° ) is the uncoupled Green's 
function of the system (probe), respectively. Indices (m,n) sums over different system eigenmodes, i.e., without the 
probing waveguides. The perturbative series can be summed to yield the exact Green's function in the form: 



where the self-energy of the coupled system is given by: E mj „ = (VGl, g V) m , n = \ Sfc uj~u k +ie V k,n- We assume 
that the group velocity of the waveguide (in the absence of the coupling to the system) is smooth around the resonators 
frequency. In this limit, the interaction coefficient can be written as: 



Vt = V , V k ° ut = V (34) 

where Vq is the coupling coefficient at the wave vector corresponding to the resonant frequency of the ring resonators 
(i.e. Vq — V m= Q t k=u /c)- We transform the momentum sum of the self-energy expression into an integral. The real 
part of the integral corresponds to a constant energy shift which can be absorbed in the energy definition. The 

imaginary part of the self-energy corresponding to input channel is: —i^\in)(in\, which shows that the field inside 

v 2 

the input resonator will decay into the input waveguide at a rate v = where v g is the group velocity at the 
resonator frequency. Similarly, one can evaluate the output self-energy terms. The self-energy for the total system 
can be written as 



S = — i— \in){in\ — i—\out)(out\. (35) 

Note that the Green's function in Eq. (f33|) . can be also written as G _1 (cj) = uj — H sys — E, which can be represented 
as a (N x N y ) x (N x N y ) matrix and easily evaluated numerically. 

Now, having the Green's function of the coupled system, we return to the evaluation of the scattered field. We 
start by evaluating the wavefunction of the transmitted field in the same input channel (as shown in Fig.(l)b in the 
main text), which can be defined as ipt(x) = lim x ^, +oc (x\4') and the input field is \k) with k > 0. One can write the 
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transmission coefficient in terms of the T-matrix: ipt{x) = 772 2fe'=-ao e%k x Tk'k- By using the T-matrix expression 
in Eq.p2p. and the Green's function, we can formally write the wavefunction of the transmitted field. The first term 
in the T-matrix gives the trivial unperturbed wavefunction -^=e lkx . The second term can be written as: 

+ °° ik' x -y 
k — — 00 

Again by transforming the momentum sum into an integral and performing the contour integral, the pre-factor before 

v 2 

the resonator summation simply reduces to — i— = —iv. Therefore, the transmission coefficient for a given frequency 

Vg 

will become: 



t(uo) = 1 - ivG in , in {u)). (36) 
Similarly, one can obtain the other transport coefficients: 

r'(uj) = -ivG outM (u k ), 
r(u) = -ivG inM (uj k ), 
t'(u) = -ivG ou t.tn{u k ). 

IV. BACK-SCATTERING DISORDERS 

The presence of scatterers in the system leads to mixing of pseudo-spin up and spin down, acting as a "magnetic 
disorder" . If this coupling is done in a controllable fashion, we arrive at the spin coupling Hamiltonians discussed 
earlier. However, in a experimental realization undesired backscatterings (e.g., due to surface roughness) lead to a 
similar coupling between spin-up and -down components. In most cases, this coupling is negligible. This has been 
the case for coupled resonators system 3^, 45-48|. However, in very high-Q resonators, the mode coupling (e.g., 



due to surface roughness) is usually greater than the resonator decay rate. In this case, the resonator enhances 



of the backscattcring effect, and a mode splitting occurs [49|, |50|. We note that as long as the interest is in the 
linear regime (as in this article), very high-Q resonators are not required. Nevertheless for completeness, we study 
the backscattering effect by first characterizing it and then evaluating its effects. As we show below, the effect of 
backscattcring is more pronounced in the resonators and can be modeled by a magnetic disorder Hamiltonian. Before 
showing the derivation of the magnetic disorder Hamiltonian, we present the result of such model. In particular, we 
consider a case where at the resonator (x! , j/), a scatterer is present and characterized by transmission (reflection) 
coefficient as t s ~ 1 (r s = respectively. The corresponding perturbation will be: 

*~=^.,.u>(A e ;)fe:) 

where <p depends on the position of the scatterer inside the resonator. The effect of such term can be understood 
in the following way: If we consider a forward-going spin-up wave, the disorder (to the first order in e) scatters the 
photon current into forward- and backward-going spin-down wave. Similar effects occur for other spin currents. This 
is in direct analogy to a "magnetic disorder" in the context of electronic QSHE, where electronic spin currents are 
not immune to magnetic impurities. Since our system is bosonic such impurities does not break the time reversal 
symmetry. However, their presence has a similar effect to magnetic impurities in fermionic QSHE. In particular, if we 
excite the system in spin-up modes, due to backscattering, spin-down modes will be also populated. However, since 
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Figure 7: Scattering of a traveling spin-up photons from a resonator with a disorder 

spin-up and -down are coupled to different output channels (i.e., (r',t) and (r, if), respectively), the backscattering 
effect will be manifested in the reduction of signal in r' , t channels and some leakage in r, t' channels. Fig. 4 in the main 
text shows different transport coefficients where we have assumed the backscattering is as strong as the dissipative loss. 
We observe that although the transport properties of magnetic bands states are affected by such magnetic impurities, 
edge state transport is robust, which is due the suppression of backscattering events. In particular, the scattering of 
a forward-going spin-up into forward (backward) -going spin-down is inhibited due to energy (momentum) mismatch. 

Here we show that the effect of having an scattcrer in a resonator can be modeled by a simple perturbation in our 
Hamiltonian formalism. The scatterer inside the resonator back-scatters the photon in the opposite mode, which in 
our spin language means that the spin-up photons can be scattered into spin down photons, vice verse. In particular, 
we consider the following scenario: a spin-up forward propagating photon current, in the presence of a back-scatterer 
in a resonator on its way, can get back-scattered into forward- and backward-going spin-down current (Fig{7j) . We 
first study this scenario in the physical system of a ID array of resonators by transfer matrix formalism to evaluate 
the forward- and backward-going spin-down field amplitudes. Then, we study a model Hamiltonian and evaluate the 
scattered field amplitudes. By comparing the physical system and the model Hamiltonian, we argue that the physical 
system can be faithfully modeled by our proposed Hamiltonian. 

First, we consider a ID array of resonators and locate a disorder inside a single resonator and evaluate the transport 
properties of the a traveling spin wave off that disorder. In the absence of disorder the system possesses four modes, 
corresponding to clockwise and counter-clockwise modes of the resonators where each can propagate forwards or 
backwards. Therefore, at each site the photonic field can be represented by a 4 dimensional vector. We work in the 
basis of this four modes: so for example (1, 0, 0, 0) T represents the forward spin-up photons and similarly for the three 
other modes. 

In this basis, as discussed in the Sec. I of this SI, the transfer matrix of each cell is diagonalized in the following 
form: 



M = - 



1 e lKA 








\ 





e -iKA 














e lKA 





V o 








e~ lKA J 



(38) 



where cos(ifA) = — ^ and < l. Now, we consider that there is a cell with a disorder located at site zero 
(...M MoMoM sca uM M ...). The propagation of a free travelling spin wave can be deduced by application of the 
transfer matrix (MoM scat tMo). In particular, we are interested in the reflected (££~) and the transmitted field £^+i- 
However, these output field are not readily available in the transfer matrix, and we have to use the S-matrix. More 



precisely, n + 
V & n+i 




while for S-matrix: 




5' 



In our case, we deal with two 



n+l 



forward- and two backward-modes, so the M-matrix is characterized by 2x2 blocks as: M 

.. o i c ( A-BD- l C BD- 1 \ 

corresponding S-matrix will be given by S — \ , , 

I -D~ X C D- 1 J 




therefore, the 
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We consider a situation where the disorder is in the middle of the resonator in the upper arm. Therefore, the 
transfer matrix of two cells with a disorder in between will be M — M M scatt M . The corresponding S-matrix to the 
hrst order in the disorder strength (e) and the zeroth order in the inverse of the finesse (J 7 ^ 1 = (1 — r 2 )/-7r), will be: 



2iKA 



( 1 o o\ 
10 
10 

y o o o i j 



2iKA 



( 1 1 \ 
10 10 
10 1 

V i o i o j 



(39) 



where the first part simply represents the free propagation over two cells, and the second part represents the scattering 
into the other modes due backscattering. Note that as the scatterer is in the resonator, its strength is enhanced by 
the finesses of the resonator. Now, we consider a travelling spin-up wave which goes forward. For convenience, we 
assume the phase of this wave is such that it is zero at the site with the disorder. Therefore, our input will be: 



-iKA 






\ 



. Since we are considering the propagation over two cells, in the absence of the scatterer, the 



V o J 

output field should be: v € ^ ut = S(e — > 0)vi. 



( e+* KA \ 



V o / 



However, in the presence of the scatterer, we will have: 



Vout 



,+tKA 



1 \ 

7r sin(KA) 

o 



\ Trsm(KA) I 



The second (forth) row terms correspond to spin-down forward- (backward-) going 



amplitudes. 

Now, we consider a tight-binding model with a disorder and show that it has the same behavior as our photonic 
system. In particular, we assume that the Hamiltonian of the system to be: 



H 



-E(«i«i)(;;)(:)-^4 




(40) 



where the first term corresponds to a simple spin-conserving hopping Hamiltonian and the second terms corresponds 
to a magnetic disorder at site zero which rotates the pseudo-spin along the x-axis at the specific site. We assume the 
disorder to be weak (e' <C 1). Since the system is linear, the dynamics of the system can be studied by considering 
a single excitation. Therefore, a general state of the system can be captured in a two-component wavefunction as a 
function of sites: |\E f ) = J2 m ^f,i( m ) a | JO). In particular, since we are interested in the scattering of a forward-going 
spin-up excitation off a disorder positioned at site m=0, the wavefunction of the entire system can be described as: 



e iKAm _|_ -iJCAm 



1^(0) 

+iK Am 



m< 

m = 
m > 



(41) 
(42) 
(43) 



ne -iKAm m < 

^4.(0) m = 

I t ie +iKAm to > 



(44) 
(45) 
(46) 
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where 4-f (tj.), corresponds to the transmitted amplitude of the spin- up (-down) component. Similarly, r-f (r±), 
corresponds to the reflected amplitude of the spin-up (-down) component. Now, using the above ansatz, we solve the 
Schrodinger equation: H\^ K ) — E K \^ K ), for m = 0,±1 and \m\ > 1, which leads to seven independent equations: 



E k = -2Jcos(KA) (47) 

E k ^(0) = -J(e- lKA + r t e* KA )-Jt t e +lKA -e'JiP i (0) (48) 

EMO) = - J n e iKA -Jt i e+ tKA -e'J^(0) (49) 

E k t t e lKA = -J^ t (0) - Jt f e +t2KA (50) 

E k t ie lKA = -J^(0)- Jt i e +aKA (51) 

E k (e- tKA + r t e iKA ) = -J^(0)-J{e~ 2iKA + r t e+ t2KA ) (52) 

E kn e lKA = -J^{0)-Jne +l2KA . (53) 



We can solve the above equations for an arbitrary value of e' . However, since we are only interested in weak impurities, 
we present the solutions to the first order in e': 



t t = 1 , r r = (54) 

^ = 1 2^ka~y (55) 

These results mean that a forward-going spin-up wave scatters into forward- and backward-going spin-down wave. 
This is identical to the perturbative result of our physical system (i.e. resonators and waveguides). In particular, if we 
set J = k and e' = — 2eF/ii, then the model Hamiltonian in Eg I3TI and Eq. 001 are equivalent to each other and they 
describe the scattering inside the resonator. We note the the above derivations can be also reproduced for scattering 
of a spin-down wave propagating in the forward or backward direction. 

Moreover, we note that where the scatterer is not positioned in the middle of the resonator the effective Hamiltonian 
is no longer an in-plane magnetic field in the x-direction. For a general position, the model Hamiltonian takes the 
form: 



— e«4,)(;: )(:;;) -v(4o %•*)(::) 

where the scattering terms corresponds to a magnetic field pointing in a random direction in the system plane. 

Furthermore, we note that one might be tempted to try an educated guess based on the spin Hall Hamiltonian 
presented in SecJITl without performing the above calculations. However, due the high symmetry of the system 
considered in Sec|TT] (symmetric positioning of the scatterers), we can not arrive at a rigorous approach to find the 
exact form and the coefficients of Eq. (|4"U)) . In particular, varying the position of the scatterer in the single disorder case 
leads to a magnetic field pointing in any direction in the plane, while in the spin Hall model of SecHJ symmetrically 
varying the position of the scatterers leads to only a strength modification of the in-plane magnetic field. 



EFFECT OF LOSS 



Here we discuss the effect of photonic loss which leads to leakage of photons out of the system. The loss can occur 
both in resonators and in the connecting waveguides. In most system, propagation and bending losses dominate the 
light attenuation and one can ignore losses due to coupling regions 33|, l5l| . We characterize the propagation loss and 
the bending loss at 90° turns by a single quantity /i which is the field attenuation in a roundtrip inside the resonator. 
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By assuming moderate quality factor for resonators and using transfer matrix formalism, one can easily show that 
the dispersion relation of a ID array of resonators takes the form: 

2 

ui(K) = ojq — 2k cos(A'A + 2noiy) — i—^Kj 7 . (57) 

7T 

This loss can be also represented by a non-hermitian term in the Hamiltonian: 

Hi oss = -iK in a\a x (58) 

where the dissipation rate is given as the product of loss per roundtrip, number of roundtrips per decay time and the 
photon decay rate (n in — ju(^)(4k)). 

Fig. 4(c) of the main text shows the reflection spectrum in the presence of such loss. In the edge state transfer band, 
the reflection is slightly attenuated due to the attenuation of the light field during the propagation around the system 
perimeter (~ exp(— AN x Ki n )). We note that the numerical results corroborate with this estimate. 



VI. EXPERIMENTAL CONSIDERATIONS 



We next consider experimental issues involving the realization and detection of the magnetic photonic states. 
Although our formalism is general and can be applied to various photonic systems, here, we focus on silicon-on- 
insulator (SOI) technology where it has been shown that up to 100 micro-ring resonators can be successfully coupled 
to each other operating at 1.5/zm 33[. We consider a 10x10 lattice and take the unit cell as 20/im. Since for a typical 
ring of a few microns, bending losses (~ 0.004dB/turn) are typically larger than propagation losses over a wavelength 
(~ 3.5dB/cm) 33, |5l|, using waveguides with different lengths to connect the resonators might be experimentally 



more accessible than coupling rings as in a Mach-Zender configuration [24|, [25J. Using the experimental parameter 
in Ref. jjj^], we find that the extrinsic decay is k = lnm(j^), which is half of the bandwidth of the tight-binding 
model. The intrinsic loss is n in ~ 0.02nm(j a ) = 0.02k and the standard deviation of individual resonances is 
U ~ 0.4nm(j a ) = 0.4k. In numerical simulations for Fig. 4 and 5 of the main text, we use these parameters and 
assume that the backscattering strength is e T = 0.1. 



